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We study nuclear spin dynamics in a quantum dot close to the conditions of electron spin reso- 
nance. We show that at small frequency mismatch the nuclear field detunes the resonance. Remark- 
ably, at larger frequency mismatch its effect is opposite: The nuclear system is bistable, and in one 
of the stable states the field accurately tunes the electron spin splitting to resonance. In this state 
the nuclear field fluctuations are strongly suppressed and nuclear spin relaxation is accelerated. 



Electrons confined in semiconductor quantum dots are 
being investigated intensively in recent years. Much re- 
search is inspired by the possibility to use their spin 
to implement qubits^i.e. the computational units in a 
quantum computer [lj. At present, the main obstacle 
for this development is the short spin coherence time T 2 * 
in high-purity quantum dots, measured to be in the ns 
range 0, 0] • Hyperfme coupling of the electron spin to 
randomly fluctuating nuclear spins was identified to be 
the main source of this fast decoherence [2, |4| • 

It was shown recently that hyperfme interaction in 
semiconductor quantum dots can lead to much richer 
physics than just dephasing. Various experiments have 
demonstrated a set of novel phenomena: random switch- 
ing between two stable states Q, current oscillations on 
a time scale of minutes §| , and strong hysteresis @ . All 
new effects were attributed to hyperfme induced dynam- 
ical nuclear spin polarization (DNSP) resulting from a 
non-equilibrium electron spin polarization. The nuclear 
polarization built up then feeds back to the electron spin 
splitting and is thereby observed. Optical excitation of 
quantum dots exhibited fine mode locking at multiples 
of the electron spin resonance frequency [7| ■ The large 
magnitude of the signal and slow dynamics suggest that 
DNSP tunes the resonance in individual dots [8(. Nu- 
merical simulations seem to support this point of view, 
but do not immediately supply a comprehensive picture 
of the underlying physics [9(] . 

These findings tri gger ed ideas to try to make use of 
the strong feedback |10j, and several ex per iments were 
designed to optimize the effect of DNSP [T3 | . It was 
observed that in a polarized state the fluctuations of the 
nuclear field are strongly suppressed and their relaxation 
is accelerated. Both effects may result in a significant im- 
provement of T 2 *- On the other hand, controlled DNSP 
might open up the possibility of processing quantum in- 
formation in robust nuclear spin ensembles [1 31 ] . All this 
feeds intensive research on the coupled electron-nuclear 
spin dynamics in quantum dots 
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Recently, electron spin resonance (ESR) in a double 
quantum dot [HI], has been used to perform single elec- 
tron spin rotations. Besides the demonstration of spin ro- 
tations, the experiment gave a clear indication of DNSP. 



The resonant response extended to a rather broad fre- 
quency interval without any amplitude reduction. Re- 
markably, this broadening was asymmetric with respect 
to the ESR frequency gfiBBo/h. Besides, a significant 
hysteresis was observed: The response depended on the 
sweep direction of frequency or magnetic field, suggesting 
that the resonant condition is shifted during the sweep, 
as if something tunes the electron splitting. While these 
effects have been speculatively attributed to DNSP, their 
exact mechanisms are not yet understood. 

In this work we present a simple model to study the 
coupled electron-nuclear spin dynamics in a single quan- 
tum dot close to the ESR condition, assuming that the 
ac driving is sufficiently strong to saturate the resonance. 
We find that ESR polarizes the nuclei in a preferred di- 
rection. We show that this results in tuning as well as 
detuning of the resonance by DNSP, depending on the 
mismatch between the driving frequency and gfig\Bo\/h: 
At small mismatch, the nuclear field built up simply de- 
tunes the spin splitting away from ESR. At larger mis- 
match, competition between ESR pumping and nuclear 
spin relaxation causes a bistability, and in one of the sta- 
ble states the nuclear field actually tunes the system back 
to ESR. With this model, we find a number of recently 
observed effects (strong asymmetric hysteresis @, and 
a reduction of nuclear field fluctuations and accelerated 
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FIG. 1: A quantum dot under ESR conditions. Left: The 
spin-split levels |+) and | — ) in the rotating wave approxi- 
mation. The arrows show the first-order transitions with a 
dissipative electron spin-flip (thick) and second-order transi- 
tions with a nuclear spin-flip (thin). Right: The initial and 
final states of the transitions for the case when dissipation is 
dominated by co-tunneling to the leads. 
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dynamics [l2j]), and we provide a clear explanation of the 
mechanisms involved. To achieve a quantitative agree- 
ment with the experiments, one would have to use more 
detailed and specific models to account for, e.g. the pres- 
ence of two coupled dots [15(] or the possibility of sub- 
stantial electric contributions to the ESR signal (lit - 

In our model, a single electron is confined in a quantum 
dot with its energy well below the Fermi levels of the 
nearby leads. To achieve ESR, one combines a dc and 
ac magnetic field, B to t = B z + £?i cos(wi)i, the ESR 
condition being Huj = gpB\Bo\. The interaction of the 
total magnetic field and the electron spin S is given by 
the Hamiltonian H = — g/isB to t ■ S. 

It is natural to assume that the ESR frequency mis- 
match Af = <jj — gnB\Bo\/% as well as the Rabi frequency 
B\ = gpsBi/h are much smaller than lu. This justifies 
a rotating wave approximation. In a rotating frame, we 
can write the Hamiltonian as 



H = h(Af)S z + ^BiS x 



(1) 



This Hq determines the effective electron spin eigen- 
states in the rotating frame, |+) = cos ^9\ |) + sin ^9\ [) 
and |-) = sin±0| |) - cos§0| |), where 9 = tt/2 + 
arctan(2A//Bi). To determine the probabilities p± to 
be in one of the eigenstates, we have to take into account 
dissipative processes accompanied by a spin- flip. Those 
are due to the coupling to the environment, that very 
generally can be represented as 
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H~S + e 1 ' 



H+S~e- 
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where H z,± represent the environmental degrees of free- 
dom coupled to corresponding electron spin components. 
It is convenient to consider the environment in the 
rest frame, while S is defined in the rotating frame. 
This makes the coupling explicitly time-dependent. For 
quantum dots, relevant dissipation mechanisms are: (i) 
electron-hole pair creation in a nearby lead by a co- 
tunneling process [I?}) (ii) spin-orbit interaction with 
phonons [lq |. (iii) direct coupling to phonons perturbed 
by Zeeman energy. Mechanisms (i) and (ii) can lead to 
an electron spin-flip in the z-direction in first order, not 
causing any nuclear spin flips. Mechanism (iii) couples 
only to S z and therefore it can flip the electron spin only 
through a second-order process involving the nuclei [l9j . 

We assume temperatures lower than fiui. We then find 
from energy consideration that the terms proportional 
to H~S + e luJt dominate the transition rates between |+) 
and |— ), since they correspond to the maximum energy 
transfer ps gps\BQ\ ~ Hu> from the dot to the envi- 
ronment (see Fig. [1]). To find the steady state proba- 
bilities p+ and p_, we compute these rates and solve 
the master equation T~[w] (— M_ + p + + M + _p_) = 0, 
where T~ is the maximum relaxation rate, T~ \u>] — 



${H-(0)H+{T))e~ lUT dT/ATi 2 . For mechanism (i) T~ cx 
w, while for (ii) T~ cx ui 5 . The matrix elements M a p = 
| (alls'" 1 " |/3) | 2 are calculated from |+) and |— ), and this 
yields p± — \ ± cos 9/(1 + cos 2 6). Far from the reso- 
nance (if 9 — ► 0, 7r) the spin is in the ground state 1 1) 
(corresponding with p + — ► 1 or p_ — * 1, see Fig. [T]), 
while exactly at resonance (9 = 7r/2) one finds p± = 1/2. 
This approach is valid provided that the energy splitting 
in the rotating frame is sufficiently big, i.e. B\ 3> r~[w]. 

Let us now consider the nuclear spins I, which are cou- 
pled to the electron spin via hyperfine interaction [2o| . 
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where the sum runs over all nuclei. For simplicity we 
assume that all nuclear spins are equally strongly coupled 
to the electron spin, so that the prefactor reduces to the 
hyperfine coupling energy (for bulk GaAs IE n ~ 100 
peV [i(|) divided by the effective number of nuclei N. 
For quantum dots, this number is big (JV ~ 10 6 ), and this 
distinguishes the situation in quantum dots from that of 
a single paramagnetic ion [2lj |. 

The effect of the hyperfine interaction is twofold. 
Firstly, the nuclei affect the electron dynamics via the 
Overhauser field (I z )E n , that adds to the static z- 
componcnt of the external magnetic field. Secondly, the 
interaction can cause electron spin flips to be accompa- 
nied by flips of nuclear spins. If there is a preferential 
direction for these flips, they can be the source of DNSP. 

Let us evaluate the rate of the hyperfine induced nu- 
clear spin flips. We keep in mind that hyperfine inter- 
action by itself cannot cause spin exchange between the 
electron and the nuclear system owing to the energy mis- 
match ps fiui between the initial and final state. The rate 
thus arises from a second-order process, the correspond- 
ing amplitude incorporating H coup and Hm- In principle, 
there are six processes capable of flipping nuclear spins. 
To estimate their relative magnitude, we note that the en- 
vironment favors large positive energy absorption. This 
brings us to the conclusion that the dominant process 
comes from combination of S + e luJt I~ in H^f, and S Z H Z 
in H conp (see Fig.Q]). The corresponding energy transfer 
is « Hu>. This means that nuclear spins only flip from the 
"up" to the "down" state: There is a preferential direc- 
tion, needed for DNSP [13] ■ The resulting polarization is 
negative, P = (iVf —N±)/N < 0, iV-m) being the number 
of nuclei with spin "up" ("down"). 

The pumping rate is proportional to r^[w] = 
J{H§(0)H*(T))e- luJT dT/h 2 , accounting for the dissipa- 
tion of hu>. For mechanism (i), T~ = T z owing to SU(2) 
symmetry in spin space. For spin-orbit mechanism (ii), 
r,T and are of the same order of magnitude [l8| • The 
total pumping rate reads 
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We assume P <C 1, so that Af — N/2 (even small polar- 
izations are enough to (de)tune the resonance in a wide 
frequency range) . At this stage we incorporate the effect 
of the Overhauser field: It is a simple shift of the fre- 
quency mismatch, and we define the resulting mismatch 
A = A/ + IE n P/h. We note that the validity of the 
rotating wave approximation now requires |A| -C to. The 
matrix elements in <[4|) and p± are now functions of A 
and the pumping rate assumes a Lorentzian shape 
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with the same width as e.g. (5 Z (A)). The numerical fac- 
tor accounts for I = 3/2 for GaAs. We see that T p <C T r , 
provided the Zeeman splitting tko exceeds the typical 
fluctuations of the nuclear field Kl = IE n /y/N. This 
sets the limits of validity of our perturbative approach. 

The resulting nuclear polarization follows from the 
competition between spin pumping and intrinsic nuclear 
spin relaxation characterized by the rate l/r„. In terms 
of P, the balance equation reads 



dP _ 2r p (A) 1 

dt N T n ' 



(6) 



which is in combination with the Lorentzian in Eq. ([5|) 
the main result of our work. 

To proceed, let us note that the nuclear relax- 
ation rate is very low, t„ ~ 10 s 23]. We express 
this smallness in terms of a dimensionless parameter 

a = (IS/Sv^Xt^M)" 1 ^^ 2 /^ 3 )- Asa « 
Bi/Tp(0)T n , a small a means that T p is sharply peaked 
compared to the slow relaxation rate l/r ra . Although the 
theory outlined is valid for any a, a strong DNSP feed- 
back requires a€ 1. We assume this from now on. 

Let us consider detuning first. The natural measure for 
the frequency mismatch is the resonance width B x /^2. If 
the initial frequency mismatch is small, Af / B\ < 1, then 
the weak relaxation stops the building of nuclear polar- 
ization only at significantly large A/ B\ — 2~ 3 / 2 a~ 1 t 3 3> 
1, that is, far from the resonance (see Fig. 2a). 

Counterintuitively, a larger frequency mismatch results 
in tuning. If Af/Bx > 3 • 2~ 13 / 6 a- 1/3 , then Eq. © has 
three zeros, and the polarization becomes bistable, as in 
Fig. 2b. This bistability is preserved till much bigger mis- 
matches, with upper boundary (A/) max = _Bi2~ 3 / 2 a _1 . 
In one of the stable configurations A ~ B\ , and the sys- 
tem is tuned to resonance. The other stable state is un- 
polarized and thus far away from the ESR condition. We 
stress that the bistability is asymmetric: If E n > 0, as 
in GaAs (an antiparallel arrangement of electron and nu- 
clear spins is energetically favorable), the bistability oc- 
curs only at A/ > 0. If on the other hand E n < 0, it 
occurs only at to lower than gp,B\Bo\/h. 

Such a bistability implies also hysteretic behavior. Let 
us adiabatically sweep the frequency starting at to < 
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FIG. 2: Two plots of dP/dt normalized by the maximum ESR 
pumping rate P m ax = 2r p (0)/AT, versus the nuclear polariza- 
tion P. The circles indicate the stable configurations, (a) 
Detuning at small frequency mismatch (A/ = 0). (b) Bista- 
bility and tuning at large frequency mismatch (A/ = 100-Bi). 
For both plots a = 10" 3 . 



g^B\B{)\ far from the resonance (see Fig. 3). Upon in- 
creasing to we first cross <?/xb|-Bo| (i.e. the line A/ = 0), 
and then get to the tuned state at A/ ~ B\oT x l 3 . 
We remain in this state until the frequency mismatch 
reaches (A/) max , provided that our sweep speed is much 
smaller than the typical nuclear spin pumping rate, i.e. 
to <C 2E n T p (0)/N'h. If we cross (A/) max , nuclear relax- 
ation becomes stronger than DNSP: The tuning ceases 
and the only stable state is again the unpolarized one, 
which will be reached on a time scale of t„ . If we then go 
backwards decreasing to, we will not get into the tuned 
state but remain in the stable unpolarized state. 

The overall structure in the (Bq, o;)-plane is sketched in 
Fig. 3, where the bistability occurs in the gray-shaded re- 
gion. An experimentally accessible quantity is the width 
of this region. It is equal to the maximum frequency 
mismatch (A/) max oc a -1 a r^[w]/w 2 , so it exhibits 
dependence on to. At lower frequencies, mechanism (i) 
(interaction with the leads) dominates the dissipative 
spin-flips and r*[o;] oc to, so that (A/) max oc 1/to. At 
larger frequencies mechanism (ii) takes over, resulting 
in (A/) max oc to 3 . This together implies that (A/) max 
reaches a minimum at a finite to c . We illustrate this be- 
havior in Fig. 3, choosing (A/) max = 0.04o; c [3(w c /w) + 
(lo/uj c ) 3 ]. The separate contributions of mechanism (i) 
and (ii) are indicated with thin dashed lines. 

Let us give a numerical example supporting the as- 
sumptions made. Based on typical experimental param- 
eters [H, we use T n = 15 s, N = 10 6 , B x = 1.5 mT, 
to = 120 mT and Q = 5 mT (for \g\ = 0.35 15], 1 mT 
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We take Tl = 2 ■ 10 b s"\ this is in ac- 



cordance with a lower bound estimate set by the typical 
leakage current of 100 fA [lif. We find that a « 1.5-10" 2 , 
so it is small indeed, and this suggests strong DNSP. The 
same set of parameters gives (A/) max w 24Bi, which is 
much bigger than the resonance width and even compa- 
rable to the driving frequency (A/) max rj 0.3 a;. This 
gives 40 mT, a value that agrees well with experimental 
observations in a double dot ([IB], Fig. 2C). The maxi- 
mum polarization is achieved at the edge of the region, 
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where f(q) = arctan \/ q~ x — 1 — ^/qJT~~q) ~ 1. The 
values from our numerical example give 1 /r t ^ u ~ 4 s at 
g = i. The inverse rate r„_^ t is yet smaller, so that if the 
dot has switched to the unpolarized state, it is unlikely 
to switch back by itself. One would have to make again 
a frequency sweep as described above. 



FIG. 3: Stability diagram in the (w, Bo)-plane for ESR in- 
duced DNSP. Bistability and tuning to the resonance occur 
in the gray-shaded region. The arrows show an adiabatic fre- 
quency sweep leading to the tuned state. The dashed lines 
present the separate contributions of mechanism (i) and (ii). 
This diagram qualitatively agrees with [TJj], Fig. 2C 



where P max w 7.2 • 10~ 3 , so that P < 1 indeed. 

We also investigated both the switching rates between 
the tuned and the unpolarized state, and the small fluc- 
tuations near these stable states. To estimate the fluc- 
tuations, we use a Fokker-Planck equation for the distri- 
bution function of the polarization P(P). To derive the 
equation, we regard the nuclear dynamics as a random 
walk on a discrete set of spin values n=l(N< [ -N i ). The 
pumping rate T p only causes transitions from n to n — I, 
while the spin relaxation causes transitions in both di- 
rections with almost equal rates (I/2r„)7V|,| ^ T p . We 
go to the continuous limit, justified by the large number 
of nuclei (~ 10 6 ) to obtain [24 1 
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From the steady state solution of ([7]) we evaluate the 
small fluctuations around the unpolarized and tuned 
states. While ((AP) 2 ) unp = 1/N is not affected by ESR, 
the fluctuations in the tuned state are suppressed roughly 
by a factor a -1 (67 for the numerical example), more 
precisely <(AP) 2 ) tun = \a{q\l - q)}^ ((AP) 2 ) unp , 
where q = A//(A/) max ~ 1. Importantly, this factor 
also determines the acceleration of the nuclear dynamics: 
The local nuclear spin relaxation time in the tuned state, 
T* un , is shorter than r„ by the same factor. 

This quenching of the fluctuations is also a justifi- 
cation for neglecting them. If the fluctuations would 
have been fully developed, one could only neglect them 
if the resonance width P>\ ^> f2. Since the fluctuations 
are suppressed, this condition is now achieved at much 
smaller driving fields B\ ^> fly/a. The same condi- 
tion guarantees that spontaneous switching between the 
tuned and unpolarized state occurs with an exponentially 
small rate. We evaluate this rate from ((7]) with Kramers 
method [24f to obtain 



To conclude, we have investigated DNSP under ESR 
conditions to find both detuning and bistability-related 
tuning of the resonance. The simple model in use ex- 
plains qualitatively a set of recent experimental findings. 
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